Dissipative preparation of entanglement in optical cavities 
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We propose a novel scheme for the preparation of a maximally entangled state of two atoms in an optical 
cavity. Starting from an arbitrary initial state, a singlet state is prepared as the unique fixed point of a dissipative 
quantum dynamical process. In our scheme, cavity decay is no longer undesirable, but plays an integral part 
in the dynamics. As a result, we get a qualitative improvement in the scaling of the fidelity with the cavity 
parameters. Our analysis indicates that dissipative state preparation is more than just a new conceptual approach, 
but can allow for significant improvement as compared to preparation protocols based on coherent unitary 
dynamics. 
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Entanglement has for some time now been identified as the 
key resource in quantum information tasks. In particular, for 
bipartite systems the properties and uses of maximally entan- 
gled states (Bell states) are well understood Preparing 
these Bell states faithfully and reliably has been one of the 
major challenges in the field of experimental quantum infor- 
mation science, where a plethora of different systems has been 
investigated [2]. In particular, several schemes based on cavity 
QED have been proposed (see e.g. [3-9]), and these schemes 
have been used to generate entanglement of atoms using mi- 
crowave cavities lflo[[Tlll . 

Traditionally, it has been assumed that noise can only have 
detrimental effects in quantum information processing. Re- 
cently, however, it has been suggested Il2l4l5ll . and realized 
experimentally that the environment can be used as a 
resource. In particular, it was shown in Ref. [12] that univer- 
sal quantum computation is possible using only dissipation, 
and that a very large class of states, known as Tensor Prod- 




uct States [ 171 11811 . can be prepared efficiently. On general 
grounds, one may argue that dissipative state preparation can 
have significant advantages over other state preparation meth- 
ods, but whether it really leads to an advantage for a given 
task can only be determined by considering concrete physical 
examples. Below, we argue that this is indeed the case for the 
system we study. 

In this letter, we suggest a dissipative scheme for preparing 
a maximally entangled state of two A-atoms in a high finesse 
optical cavity 1 19142 ill , with detunings as depicted in Fig. la. 
Our scheme can be understood from Fig. lb, which describes 
the effective coupled ground states of the atoms in the cavity. 
A microwave field shuffles the three triplet state around while 
an effective jump operator describing cavity leakage causes 
the transitions: |00> \S) -> where \S) = (|01) - 
\10))/y/2 is the singlet state, and |0), |1) denote logical qubit 
states. Setting the cavity detuning equal to the single atom 
cavity line shift will greatly favor the transition to the singlet 
state and strongly suppress the transition away from it, as the 
latter is far off resonance; thus driving essentially all of the 
population into the singlet. 

Our protocol actively exploits the cavity decay to drive the 



FIG. 1 : (Color online) (a) Level diagram of a single atom with laser 
detuning A and cavity detuning S from two photon resonance. The 
optical pumping laser for the two atoms differs by a relative phase 
of 7T. (b) The effective two qubit system. The driving Qmw causes 
rapid transitions between the three triplet states. The decay s e ff takes 
1 00) to \S) and \S) to |11), but as each atom in state |1) creates 
a shift in the cavity frequency by an amount g 2 /A, choosing 5 = 
g 2 /A 3> n ensures that only the [00) — > \S) transition is resonant, 
because the \S) — > |11) transition is shifted by twice that amount. 
This guarantees that | S) is the only stationary state of the dynamical 
system up to very high fidelity. The spontaneous emission rates 7 e ff,i 
will tend to reduce the fidelity by redistributing information to the 
triplet states, (c) Fidelity as a function of the cooperativity C = 
g 2 /kj. The inset gives a more accurate account of the scaling (1 — 
F < 3.5C - ) for different values of n/'y. 



system to a maximally entangled stationary state. The only 
generic source of noise left in the system is then the one com- 
ing from the spontaneous emission. This leads, quite remark- 
ably, to a linear scaling of the fidelity with the cooperativity 
(see the inset in Fig. lc), which is in contrast to schemes based 
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on controlled unitary dynamics, where there are two malevo- 
lent noise sources, cavity and atomic decay, typically resulting 
in a weaker square root scaling of the fidelity J3rfl[^]. In ad- 
dition to the asymptotic fidelity, one needs to pay particular 
attention to the time it takes to reach equilibrium. 12211 . For 
our scheme, we show that the convergence rate is rapid. 

We point out that a similar study to ours has been con- 
ducted by Wang and Schirmer I23I1 . where they consider a 
detuning of the energy levels in order to break the symmetry 
in the system, and guarantee a unique steady state. It can be 
shown that their scheme, when generalized to optical cavities, 
does not give a linear scaling of the fidelity [24], but rather 
the square root, as for coherent unitary protocols. 

The basic idea behind dissipative state preparation is to shift 
the burden of coherence and control from the unitary evolu- 
tion of the system to the effective noise from the environment. 
Here we want to engineer the system environment interactions 
such that the non-unitary quantum dynamical process relaxes 
to a pure steady state of interest, in our case the singlet state, 
in a short amount of time. 

In the following, the system-environment interaction will 
be assumed Markovian, and can thus be modeled by a master 
equation in Lindblad form: 



P = i[P, H]+Y, LjpLj - -(L^Ljp + P L]L 3 ), 



(1) 



where the Lj's are the so-called Lindblad operators. We 
derive a master equation for the ground states of the system, 
which has the singlet as unique stationary state. This is 
achieved in our setup, by constructing an effective master 
equation, whose main contributing terms are the Hamilto- 
nian H = ifiMw(-^+ + J-)> an d the Lindblad operator 
L K = y/K^\S) (00\, where J+ = |1)(0| ® 1 + 1 ® |1)(0|. 
It can readily be seen that the microwave field (JImw) 
drives the transitions between the three triplet states 
({|00), |11), \T) = (|01> + |10))/V2}), while the Lindblad 
operator, originating from the cavity field leakage, will drive 
the transitions from |00) to IS*). The singlet state is thus the 
unique fixed point of this system, and the relaxation rate 
can be seen to go as min{r2^[ W /K e ff, K e ff}> which is just the 
gap of the Liouvillian if one keeps to the weakly driven and 
strongly damped regime [25]. Other terms will contribute 
weakly to the dynamics of the system, and slightly perturb 
the stationary state away from | S) . 

We now show how to prepare this effective system in a re- 
alistic quantum optical setup. Our setup, shown in Fig. la, 
consists of two A-type three level atoms in a detuned cavity 
with two stable lower energy states |0) and |1), and an ex- 
cited state I e) with a large energy separation to the lower ly- 
ing states. We apply one far off-resonance optical laser, with 
detuning A, driving the f-> e transition and a microwave 
field driving the « 1 transition resonantly. The cavity mode 
couples the 1 O e transition off-resonantly, with detuning 



A — (5, where <5 is the cavity detuning from two photon res- 
onance. Furthermore, we assume a 7r phase difference in the 
optical laser between the two atoms. This phase difference is 
crucial in guaranteeing that the singlet is the unique stationary 
state of the reduced system. Note that one could instead add 
a phase directly into the coupling term by placing the atoms a 
half integer cavity wavelengths apart. 

In a rotating frame, this situation is described by the Hamil- 
tonian 



H 



H +H g + V+ + V- 

5(Ja + A(|e)i(e| + |e) 2 (e|) + 

[ 5 (| e ) 1 (l| + |e) 2 (l|)a + / l .c.] 

^(11)^01 + 11)2(01)+^. 
^(1^(01 -|e) 2 <0|), 



(2) 

(3) 
(4) 

(5) 



where V- — Vl, g is the cavity coupling constant, a is 
the cavity field operator, il represents the optical laser driv- 
ing strength, and 57mw the microwave driving strength. On 
top of the Hamiltonian dynamics, two sources of noise will 
inherently be present: spontaneous emission of the excited 
state of the atoms to the lower states with decay rates 7^; 
and cavity leakage at a rate n. We assume for convenience 
that the spontaneous emission rates are the same for decay- 
ing to the |0) and to the |1) states (i.e. 70 = 71 = 7/2). 
This translates into five Lindblad operators governing dissipa- 
tion: L K = ^Ra, LJ = v / ; 7 72|0)i(e|, L\ = v / 7 72|0) 2 (e|, 

Ll = ^\\)M,Ll=^J%\Ue\- 



If the optical pumping laser is sufficiently weak, and if 
the excited states are not initially populated, then the excited 
states of the atoms, as well as the excited cavity field modes, 
can be adiabatically eliminated. The resulting effective dy- 
namics will describe two two-level systems in a strongly dis- 
sipative environment. To second order in perturbation theory, 
the dynamics are then given by the effective operators [ 24] : 
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.^■-^nhV+ + V-(H^V + ) + I1, 



g (6) 
(7) 



is a non-Hermitian Hamilto- 



where H NH = H Q -\ £\ l]L 3 
nian describing the non-unitary dynamics of the excited states 
which we eliminate. Applying the above equations to our 
setup, and keeping only terms to lowest order in f2, the opera- 
tors in the effective Master equation can be evaluated explic- 
itly, yielding the effective Hamiltonian and principle Lindblad 
operator 
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H, 



eff 



r K 
^eff 



ifi M w(|l)i(0| + |1) 2 (0| + h.c.) + 0(^L) (8) 



( 5 2 /A-<5) 2 + ( K /2 + 7<5/2A) 

/ ff e W 2 
V (2 5 2 /A - S) 2 + {k/2 + 7 <5/2A) ; 



2 |5)(00| (9) 



|H)(5| 



where g e s = gfl/A. 

It can be seen from the coefficients in the principal Lind- 
blad operator, that when the cavity detuning is equal to the 
cavity line shift from a single atom, and when these are large 
compared to the cavity leakage (g 2 / A = 8 ^S> k + 7<5A), then 
the denominator of the first term in becomes very small, 
while the denominator in the second term remains large. The 
effective Lindblad operator originating from cavity leakage 
then becomes L e K ff w g effy /K/2{l/ (k/2 + 7<5/2A)|5)(00| + 
l/S\ll)(S\) « v ^ki|S')(00|. The effective Hamiltonian shuf- 
fles the triplet states among each other, so that the combined 
effect of the unitary and dissipative dynamics drives essen- 
tially all of the population to the singlet state. Hence, we have 
constructed an effective master equation which approximates 
the ideal situation described earlier. 

We now consider imperfections imposed by spontaneous 
emission. The four Lindblad operators describing sponta- 
neous emission will also transform into four independent ef- 
fective noise operators for the reduced system. In the regime 
discussed above, and keeping only terms which drive the pop- 
ulation out of the singlet state, the effective operators for spon- 
taneous emission are 



r 7 

^eff,i=l,2 



T 1 

^eff,i=3,4 




11)(5| 



\T)(S\, 



(10) 



where j eS = 7^ 2 /2A 2 . 

In order to evaluate the performance of the protocol we de- 
note by Pj the probability to be in state j, and consider the 
rate of entering and exiting the singlet state P$ = -Poo K eff,i — 
Ps( K eS,2 + 53i7eff,i)) where « e ff,i is the square of the first 
coefficient of L£ s , K e ff,2 is the square of the second, and 7^ 
are the square of the scalar coefficients of L7 S .. The effec- 
tive spontaneous emission will tend to modify the steady state 
population of the three triplet states, but we assume that Omw 
is large compared to 7 e ff and K e ff, so that all three triplet states 
are almost equally populated in steady state. Solving for the 
stationary state of the rate equation and plugging in the decay 
rates obtained from the effective operators, we get for P$ ~ 1 
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FIG. 2: (Color online) The main figure shows the population of the 
singlet state (thick dashed line) and of the three triplet states (full 
lines) as a function of time for a random initial state. The curves 
were plotted for C = 50, k = 7/2, fi = 50mw/2, g = 20fi, 
and A, S are such that they maximize the fidelity for small fi. In 
this parameter regime, the stationary state has a 92% fidelity with 
respect to the singlet state. The inset shows the maximal fidelity as a 
function of the gap size for C = 50 and k — 7/2. The main figure 
corresponds to the cross on the curve in the inset. 



Here F is the fidelity, which is given by the overlap of the 
stationary state of the dynamical process with the singlet state 
F = \(S\p„\S)\=P 3 . 

There is a trade-off between the second term in the numer- 
ator (A 2 k/2<7 2 ) reflecting the probability to generate a cavity 
photon by decaying out of the singlet state, and the term in 
the denominator (-jg 2 /2A) reflecting the scattering of cavity 
photons off the atoms. The first terms favors a small detuning 
A to increase the cavity line shift, whereas the second term 
favors a large detuning A to decrease the scattering. The op- 
timal fidelity is reached when the two terms in the numerator 
are similar (7 nA 2 /2g 2 ), in which case the two terms in 
the sum in the denominator are also similar (k w j2g 2 / A 2 ). 
This leads to an error scaling as 



1 - F oc C 



(12) 



where C — g 2 /ivy is commonly referred to as the coop- 
erativity. Plugging in the values in Eq. ( fTTT ) one gets a 
proportionality factor of roughly 3 in Eq. ( fT2l . By numer- 
ically extracting the fixed point of the full master equation, 
and then maximizing its fidelity with respect to the singlet for 
fixed values of C, we get that the actual constant is closer to 
3.5 (inset in Fig. lc), i.e. 1 — F m 3.5C -1 . This discrepancy 
can be attributed to the fact that we did not include all of 
the spontaneous emission terms in Eq. (fTTT ). In addition, the 
assumption that all three triplet states are equally populated is 
not exact. To support our analysis further, we also note that 
the fidelity scaling is essentially independent of the ratio k/j 
(inset in Fig. lc). 

For comparison, in a controlled unitary dynamics protocol, 
the fidelity will suffer errors coming from spontaneous emis- 
sion on the one hand, and from cavity decay on the other. De- 
creasing one of the error sources will typically increase the 
other in such a way that the optimal value of the fidelity is 
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1 - F oc 1/y/C H. Indeed, to the best of our knowledge, 
all entangled state preparation protocols based solely on con- 
trolled unitary dynamics scale at best as l/\/C This 
means that the linear scaling of the fidelity from Eq. (fT2] i. is 
a quadratic improvement as compared to any known closed 
system entanglement preparation protocol. We note, however, 
that it is possible to beat this if one exploits measurement and 
feedback |7H9J]. As mentioned previously, the reason for this 
improvement stems from the fact that cavity decay is used as 
a resource in our dissipative scheme, so that the only purely 
detrimental source of noise is the spontaneous emission. 

The above analysis has been conducted without any consid- 
eration of the speed of convergence. We now show, that the 
entangled stationary state can be reached rapidly. In Fig. 2, we 
simulate the dynamics of the full master equation for an ap- 
propriate set of parameters. Starting from an arbitrary initial 
state, the populations of triplet states undergo rapid coherent 
oscillations with an envelope decaying at a rate proportional to 
the gap [25], while the singlet state converges to its maximum 
value at the same rate. 

For a given cavity, g, n, and 7 can be considered fixed by 
experimental constraints, and the speed of convergence is pri- 
marily governed by the magnitudes of il and £1mw- The speed 
of convergence can be increased by increasing the driving 
laser strength (£1), but the latter can not be too large otherwise 
perturbation theory breaks down, and the excited cavity and 
atomic states can no longer be ignored. Furthermore, JImw 
can not be too small with respect to {K e ff, 7eff}, otherwise the 
coherent shuffling of the triplet states will not be sufficiently 
strong to keep them at equal population. The inset in Fig. 2 
shows how the maximal fidelity scales as a function of the gap 
for a specific set of cavity parameters. The curve is plotted by 
optimizing the fidelity, for given fixed values of the gap, with 
respect to {fi, IJmw}, for fixed values of {A, 5} (those which 
are optimal for small £1). There is clearly a trade-off between 
the accuracy of the dissipative state preparation protocol and 
the speed at which one reaches the stationary state, but close 
to optimal fidelity the dependence is weak. 

Present day experimentally achievable values for the coop- 
erativity are around C 30 lIil-21]. This puts our scheme 
at ~ 90% fidelity with respect to the singlet state. Fig. 2 
shows that the stationary state is reached in a time ~ 1000/g, 
which yields for g = (27r)35 MHz (19] a convergence time 
of roughly 5 [is starting from an arbitrary initial state. 

We have proposed a scheme for dissipatively preparing 
an entangled state of two trapped atoms in an optical cav- 
ity. From both analytical and numerical evidence, we give 
the scaling of the error explicitly, and show that the station- 
ary state is reached rapidly. Our results indicate that not only 
can one produce entanglement dissipatively in a simple cav- 
ity system, but, to the best of our knowledge, the scaling of 



the fidelity for such entanglement preparation is better than 
any existing coherent unitary protocol. These results are an 
indication that an approach based on dissipation can be very 
fruitful for state preparation, as one manifestly can transform 
a previously undesirable noise source into a resource. It would 
be interesting to see if one could obtain similar results in re- 
lated systems such as trapped ions and solid state based quan- 
tum devices, where dissipation traditionally plays a detrimen- 
tal role. 
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